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ABSTRACT 

Recent 7-ray observations suggest that the 7-ray mihisecond pulsar (MSP) pop- 
ulation is separated into two sub-classes with respect to the pair multiplicity. Here, 
we calculate the cosmic ray electron/positron spectra from MSPs. Based on the as- 
sumption of the equipartition in the pulsar wind region the typical energy of elec- 
trons/positrons ejected by a MSP with the pair multiplicity of order unity is 50 
TeV. In this case, we find that a large peak at 10 - 50 TeV energy range would be 
observed in the cosmic ray electron/positron spectrum. Even if the fraction of pair 
starved MSPs is 10%, the large peak would be detectable in the future observations. 
We also calculate the contribution from MSPs with high pair multiplicity to the elec- 
tron/positron spectrum. We suggest that if the multiplicity of dominant MSP popu- 
lation is ~ 10'^, electrons/positrons from them may contribute to the observed excess 
from the background electron/positron flux and positron fraction. 

Key words: stars: neutron — cosmic rays. 



1 INTRODUCTION 

The Fer mt Gamma-Ray S pace Telescope has detected 7-ray pulsed emissions from more than twenty millisecond pulsars 



(MSPs) (|Abdo et al.ll2011h . which have a rotation angular frequency f2 ~ 10 s and a stellar surface magnetic field Bs 



lif -^G. The detection of the GeV emissions from a pulsar magnetosphere means that electrons and positrons are accelerated 
to more than ~ TeV by the electric field parallel to t he magnetic field, which arises in a depleted region of the Goldreich- 
Julian (GJ) charge density i Goldreich fc Julian 19691 ). The 7-ray light curve is an important tool for probing the particle 
acceleration process in the pulsar m agnetosphere. Therefore, the 7-ray emiss ion region has been explored by comparing 
theoreti cal models su ch as polar cap ( Daughertv fc Hardinj 19961 ). outer gap I Cheng. Ho fc Ruderman 19861 ) and slot gap 
models ( Muslimov fc H arding 2004ah with the observed light curve (e.g.. Venter, Harding & Guillemot 2009; Romani & 
Wal ters 2010; Kisaka fc Kojima 2 11). 

Venter. Harding fc GuillemotI (|2009l ) fitted the pulse profiles of the Fermi detected MSPs with the geometries of 7-ray 



emission region predicted by different theoretical models. They found that the pulse profiles of six of eight MSPs could be fitted 
by the geometries of either the outer gap or the slot gap model, as was the case of canonical pulsars. They inter preted that 
copious pairs are produced in the magnetosphere of these MSPs. However, IVenter. Harding fc Guillemot l|2009l ) also found 
that the pulse profiles of remaining two MSPs show the unusual behavior in the 7-ray light curves and could not be fitted by 
the geometry of either the outer g ap or the slot gap models. Th ey proposed that these unusual light curves could be fitted by 
the pair starved polar cap model ijMuslimov fc Hardineil2004bl ). in which the multiplicity of the pairs is not high enough to 
completely screen the electric field above the polar cap, and the particles are continuousl y accelerated up to high altitu d e over 
the entire open field line region. Thus, from the model fitting of the 7-ray light curves. Venter. Harding fc Guillemot ( 20091 ) 
suggested that the 7-ray MSP population is separated into two sub-classes. 

The important fact is that radio pulsed emission is also detected from all currently detected 7-ray MSPs and remarkably 
similar to that from canonical pulsars. The pulsar radio emission is a highly coherent process because the brightness tempera- 
ture is extremely high. In the theoretical models of the radio emission mechanisms, some authors have believed the conditions 
that there are a highly relativistic primary beam with the large Lorentz factor (~ 10^) and the number density nearly equal 
to the GJ density, and the secondary electron/positron plasma with relatively small bulk streaming Lorentz factor (~ 10 - 
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10'^) and the large pair multiplicity (~ 10"^ - 10^) in the radio emitting region (e.g., Melrose 1995; Lyutikov, Blandford & 
Machabeli 1999; Gedalin, Gruman & Melrose 2002). However, the existence of pair starved MSPs suggests that the radio 
emission mechanisms should be insensitive to the particle number density down to sub-GJ number density. The pulsar radi o 
emission mechanism is still poorly understood, so that the observationally-based constraints are valuable ( Melros j 1995 ). 
Therefore, another verification for the extent of the MSP multiplicity, especially the existence of the pair starved MSPs is 
important for the pulsar radio emission mechanisms. 

Recently, HESS has discovered a new TeV source, which is located in the close vicini ty of the globular cluster Terzan 5 
jAbramowski et ahlboil ) . Several globular clusters, including Terzan 5 also emit GeV 7-ray ( Abdo et al.ll2009l : Kong. Hui fc Cheng 



2O10l: 



Abdo et al. 



2O10l: 



Tam et al 



2011), which may plausibly be due to a number of MSPs residing in these clusters 



I Harding. Usov fc Muslimov 2005 : Venter, de Jager fc Clapson 20091 ) . Thus, inverse Compton scattering by the high-energy 



parti cles ejected from MSPs are proposed for the origin of the observed TeV emission (jBednarek fc Sitarek 2007l : IVenter. de Jager fc Clapson 
l2009r) . The high-energy electron/positron spectrum ejected from MSPs would be a useful probe for the multiplicity of the MSPs. 
However, only from the TeV spectra, we cannot di stinguish two m odels tee dnarek fc Sitarekll2007 : Venter, de Jager fc Clapson 
2009h . which assume different pair multiplicities IjAbramowski et al. 2011,). Another way to investigate the electron/positron 



spectrum ejected from MSPs is its direct measurement. Since high-energy electrons/positrons can propagate only about a 
few kpc due to the energy losses by the synchrotron and the inverse Compton emission, the direct detection of the elec- 
trons/positrons ejected from MSPs in the globular clusters is unlikely. However, for the following reasons, we may detect 
those from nearby MSPs. 

MSPs have much lower spin-down luminosity than canonical pulsars. Biisching. Venter fc de Jager ( 20081 ) investigated 
the possible contribution of the nearby MSP, PSR J0437-4715 to the cosmic ray electron/positron spectrum. They concluded 
that unlike canonical pulsars such as Geminga pulsar, the contribution from a MSP to the observed electron/positron fiux is 
negligible. However, since the lifetime of MSPs is much longer than canonical pulsars (> 10^" yr), there should be much more 
nearby active MSPs. Furthermore, Kashiyama, loka fc Kawanaka (2011; hereafter KIKll) pointed out that since white dwarf 
pulsars have long lifetime and continue to inject the electrons/positrons after the nebulae stop expanding, the adiabatic energy 
losses of electrons/positrons in the pulsar wind nebula region are negligible. Also the synchrotron cooling of electrons/positrons 
is so small and the high-energy electrons/positrons can escape the nebulae without losing much energy. Although they consider 
the case of the white dwarfs, their results are also applicable to MSPs. Therefore MSPs could potentially contribute to the 
observed high en ergy cosmic ray elec trons/positrons and wil l be detectable by the next generation experiments, such as 
CALET jTorii et al. 2008a) and CTA (|CTA Consortiumlboiol ). 

In this paper, we investigate the contribution of electrons/positrons ejected from the MSPs to the observed cosmic ray 
spectrum. In section 2, we apply KIKll model to the case of MSPs. We estimate the typical energy of electrons/positrons from 
the MSPs and show that during the propagation in pulsar wind nebulae, the adiabatic losses and radiative cooling of elec- 
trons/positrons are not so large. We also describe the propagation in the interstellar medium (ISM). In section 3, we calculate 
the energy spectrum of cosmic ray electrons/positrons from the MSPs and show the possibility that the electrons/positrons 
from these MSPs are detectable for the future observations. 



2 THE MODEL 

2.1 Acceleration and cooling 

In order to estimate the energy of electrons/positrons available in the wind region and their adiabatic and radiative cooling 
in the shocked region, we adopt the model of KIKll. For a pulsar wind nebula formed by a MSP, we consider the conditions 
that the relativistic wind blasts off from the light cylinder ~ Ric — where c is the speed of light, and two shock fronts are 
formed between the supersonic pulsar wind and ISM. Since the energy from a MSP is continuously transported to the wind, 
the shock fronts are expanding until the pressure of the shocked region P^h becomes equal to that of ISM p. Although KIKll 
considered the case of white dwarf pulsars, the situations are similar to the case of MSPs because they have a long lifetime 
and the supernova shock front have already decayed. We assume that the effects of binary companion are negligible, because 
the fraction of solid angle occupied by companion is small (<1%). We also neglect radiative loss due to curvature radiation 
within light cylinder (~ 10%). From now on we set fiducial parameters of the MSP's surface magnetic field strength, angular 
frequency and radius as Bo = 10^'^G, Q, = lO'^s"^ and R — 10® cm, respectively. 

We assume the energy equipartition between particles and magnetic field, SeN — B'^/8tv, and the conservation of the 
particle number flux, Aivr^cN ~ constant, in the MSP wind region. Here, is the number density of electrons/positrons. The 
number density can be described as 

where k is the multiplicity of electrons/positrons, A^ic and Bi^ are the number density and the magnetic field at the light 
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cylinder, respectively. We assume the magnetic field configuration as pure dipole {B cx within light cylinder. Outside the 
light cylinder, we assume the conservation of the energy flux of the magnetic field, Br ~ constant. Using these assumptions, 
the typical energy of electrons/positrons Ee can be described as 

where ^/imax is the electric potential difference across the open magnetic field lines described as 

These values are similar to those in the case of white dwarf pulsars (KIKll). Note that the typical energy depends on the 
pair multiplicity. 

Next, we estimate the adiabatic and the radiative cooling of electrons/positrons in the shocked region. The outer shock 
of the pulsar wind nebula finally decays when the pressure of the shocked region Pah becomes equal to that of the ISM p. If 
the outer shock decaying time is shorter than the lifetime of MSP, the adiabatic cooling is negligible. In order to estimate the 
outer shock decaying timescale, we solve the equation of motion and the energy conservation law at the outer shock front. 
The equation of motion can be described as 

d j 4n 3 dRout\ . „2 n fA\ 

where -Rout is the radius of the outer shock front. Pah is the pressure of the shocked region and p is the density of ISM. The 
energy equation is 

where Lad is the spin-down luminosity of MSP, 

B'o^'R' o.,.n33/ Bo Yf n Yf R 



In the derivation of eq.©, we assume that in the shocked region the internal energy of particles is 3Psh/2 because the energy 
of particles is the relativistic regime. Using the typical value for the density of ISM p ~ 10~^*g cm~^, the pressure in ISM is 
p ~ lO^^'^dyn cm^^. Solving above equations, the outer shock decays at about 



(8) 



where T is the temperature of ISM. The lifetime of a MSP r can be estimated as 

-Erot 
-tJsd 

where Erot is the rotation energy described as 



M \ f R Y f n 



2 



-~5X10^«(-:— )(-^) {-£-] {-£-] yr. (10) 



where M is the mass of a MSP. For the fiducial parameters of MSP 

M \ f Bo \ f ^ \ f_R ' 
1.4MQy VlOS-^Cy VlO^s-V VlOCcm 
We find that the outer shock decays at a very early stage of the lifetime of MSP and does not expand after t > t^cc- Therefore, 
in the similar to the case of white dwarf pulsars (KIKll), the adiabatic cooling shall give minor contributions to the cooling 
process of the high-energy electrons/positrons. 

Also for the radiative cooling by the synchrotron radiation in the shocked region Pin < r < Pout, we follow the discussion 
in KIKll based on the diffusion in the shocked region. We take the Bohm limit, where the fiuctuation of the magnetic field 
is comparable to the coherent magnetic field strength. In this limit, the timescale tdis for the electron/positron trapping in 
the shocked region is given by 

ZeBd? 

iDah 2 EeC 

where Dah = crg/3 is the diffusion coefficient under the Bohm limit, d is the size of the shocked region and rg — Ee/eB is 
the Larmor radius of electron/positron with energy Ee- To estimate the diffusion timescale, we need to know the size and the 
magnetic field strength of the shocked region. Here we consider the time t > tdoc, so that the size of the shocked region is an 
order of the radius of forward shock front at t = tdcc, d ~ Rout{t = tdcc) '--^ 3 x lO^^cm. For the radius of the inner shock front 
at t > tdcc, we use the balance condition between the momentum transferred by wind and the pressure of ISM 

1^=P- (12) 
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For the fiducial parameters, Rin{t > tdec) ^ 2 x 10^^ cm. The strength of the magnetic field at the inner radius Bin can be 
estimated as -Bin ~ 2 x 10~^ G, which is almost the same as that of ISM. Then, the diffusion timescale is 

*di«-2xl0*r-^)"yr. (13) 



.50TeV 

The synchrotron energy loss of a particle with energy is described as 

= --(TTC/? — J , (14) 

at 3 Stt ynieC'' J 

where itt is the Thomson scattering cross section, /3 is the particle velocity normalized by the speed of light and nie is the 
mass of electron/positron. Then, The typical energy loss of the electrons /positrons Ase with energy Ce can be estimated as 

17~°-KtoSg) (toS^) (toS^) • (^^^ 

This means that the high-energy electrons/positrons injected into the shocked region lose roughly 30% of the energy by the 
synchrotron radiation before diffusing out into ISM. Therefore, as in the case of white dwarf pulsars (KIKll), we can conclude 
that the radiative energy loss of electrons/positrons in the pulsar wind nebula is not so large. 

The above expressions for the estimate of the energy losses are only applicable to the case that the veloci ty of a MSP 
is su bsonic in ISM. The observed velocity of MSPs is less than that of canonical pulsars in averag e sense |Hobbs et al 



2005h. However, som e MSPs have the large velocity and a few MSPs actually forms bow shock nebulae (jStappers et al.ll2003l : 



Hui fc Becke^bood l. In this case, the size of the bow shock is described as (e.g., Wilkin 1996) 



^•^"^ - {^^) ~ (ioSg) (io£^) (lo^) (low s-0 ^^^^ 

where V is the velocity of a MSP. Due to the assumption of the energy equipartition, the strength of the magnetic field is 

^bow = ( ^) ~ 50 ( j^r^^—r) mG. (17) 



The ratio of the Larmor radius of electrons/positrons to the bow shock radius is 

-^^0.5k-\ (18) 

-'•'bow 

The fact that rg/_Rbow is close to unity supports that electrons/positrons may escape from the bow shock region. Therefore, in 
the c ase of k ~ 1, high-energy electrons/positrons can escape with an efficiency of order unity I Bandiera 20081 : Blasi &: Amatol 
2O10h . Even if we consider the case of k 3> 1, the synchrotron loss can be estimated by using eqs. (11), (14) and (16) as 



Aee „ Bo y f Q Y f R Y f V 



Ee -^""^^ (v10«-"gJ (vlO^s-iJ \m^cxa) \Wcms-^)- ^^^^ 
Therefore, we can conclude again that the radiative energy loss of electrons/positrons in the pulsar wind nebula is not so 
large. 



2.2 Diffusion in Interstellar medium 



The observed electron/positron spectrum after the propagation in ISM is obtained by solving the diffusion equation 

^/(t, r, Ee) = D{Ee)S7\f + ^ {P{ee)f) + Q{t, T, Ee), (20) 

where fit, r, Ee) is the energy distribution function of electrons/positrons, D{£e) = Do(l+£e/3GeV)'' is the diffusion coefficient, 
P{Ee) is the cooling function of the electrons/positrons which takes into account synchrotron emissions and inverse Compton 
scatterings during the propagation, and Q{t,Ee,r) is the injection term. Here we adopt Do = 5.8 x 10^*cm^s~^, 5 = 1/3, 
which is consistent with the boron-to-carbon ratio according to the latest GALPROP code. Atoyan et al. (1995) showed a 
solution in the case of an instantaneous injection from a single point-like source, i.e. Q{t,Ee,r) ~ Qo{Ee)S{t — ti)5{r). Then 
the observed spectrum G{t,r,£e;i) would be 



G{t,r,Ee;i) 



Qo(ee,o)f (ee,o) 



■ exp 



(21) 



7r3/2p(eg)ddiff(ee,ee,o)^ "'^ V ddiB{£e,Ee,o)^ 

where Ee,o is the energy of electrons/positrons at the time t{< t) and which are cooled down to £e at the time t, and ddiff is 
the diffusion length given by 



ddiff = 2 



The cooling function P{Ee) can be described as 



AaTEi 



Stt 



P{x) 

dEjUtot{Ej)fKN 



(22) 



(23) 
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where Utot{£-t)de~, is the energy density of interstellar radiation fields with the photon energy between and e~f + de~f 
(including cosmic microwave background, starlight, and dust emission; Porter et al. 2008), and B is the interstellar magnetic 
field with we here set as fpG. Here the function fK^{x) is the correction factor to include the Klein-Nishina (KN) effect, 
which approaches to unity when x is much smaller than unity. The mathematical expression of /kn can be found in Moderski 
et al. (2005). 

As shown in the last section, in general MSPs have a very long lifetime (r ~ 5 x 10^°yr) which is comparable with the 
cosmic age. In such a case that a point-like source with a finite duration, taking into account the time-dependence of an 
injection rate (Qo(£e) Qo{£e,i)), the spectrum can be calculated by integrating G{t,r,ee;T) for t: 



fi{t,r,ee;ti) ^ di G {t,r,ee;t) , 



(24) 



Qo(£e, i) oc e-'^ exp ( ) ( 1 + ^ ) , (26) 



where ti is the time when the particle injection from a source has started. Here we assume that the electron/positron injection 
spectrum can be described as mono-energetic distribution 

Qo(£e,t)oc M + ^j , (25) 

or power-law distribution 

Qo(£e, i) (X e~°' exp ( - 

Ecut ^ 

where a is the intrinsic power-law index of an electron/positron spectrum and ecut is the maximum electron/positron energy 
from a source. 

Now we can calculate the average electron/positron spectrum by considering the birth rate of MSPs as follows: 

/avo(£e) = / dti / 2Tvrdrfi{to,r,£e:,ti)R, (27) 

Jo Jo 

where to is the cosmic age (i.e. the present time), and R is the local pulsar birth rate (yr~^ kpc~^). Here £e,i is the energy at 
the time U of CR electrons/positrons which are cooled down to Se at t. 



3 RESULTS AND DISCUSSIONS 

First, we calculate the cosmic ray electron/positron spectra from the pair starved MSPs. We set the pair multiplicity k = 1, the 
lifetime r = 5 x 10^"yr, the total energy Erot ~ W^^ erg, the local birth rate R = 3 x 10~^ yr~^ kpc~^ and the fraction of the 
lost energy due to synchrotron emission 30% for each MSP. We assume that each MSP has the same value of the parameters 
(Bo = 10*-^G,n = 10^s'^,R = lO^'cm), because most MSPs have the almost same spin-down luminosity. For the injection 
distribution function, we assume mono-energetic distribution ea. (|25|) with the energy Se = 50 TeV. Even if we consider the 
power-law distribution ea. (|26p . the energy range of the distribution is small because the allowed ra nge of the c utoff energy 



should be only 50 - 80 TeV due to the observed constraint by KASKADE/GRAPES/CASA-MIA (jKistler fc Yuksel 2009 
Thus the distribution should be nearly mono-energetic distribution. Note that our local MS P birth rate is based on the MSP 
local surface density of 38± 16 pulsars kpc~^ for 430 MHz luminosity above 1 mjy kpc^ (|Lorimer"2008'). Actually, now 23 



MSPs are detected within Ikpc ( ATNF catalog. Manchester et al. 20051 '). We can only detect MSPs that have the radio beam 



directed toward us and the radio flux larger than the threshold of detectors. However, the cosmic ray electrons/positrons 
ejected from MSPs are distributed isotropically because of the effect of Galactic magnetic field, so that a large number of 
MSPs will contribute to the observed electron/positron spectrum. Therefore, our local MSP birth rate corresponds to the 
lower limit for the current radio observations. 

In figure [T] the electron/positron flux from multiple pair starved MSPs is shown. Thick solid, dashed and dash-dotted 
lines show the total electron/positron spectra if the fraction of the pair starved MSPs is 100%, 25% and 10%, respectively. 
Thin solid line shows the contribution of electrons/positrons from the pair starved MSPs when the fraction is 100%. The 
background flux is shown as a thin dashed line. We adopt the background mo del of an exponen tially cutoff power law with 



an index of -3.0 and a cutoff at 1.5 TeV, which is similar to that shown in lAharonian et al.l (2008) and reproduces the 
data in ~10 GeV-1 TeV well. It is very interesting that there is a large peak at 10-50 TeV energy range. The existence of 
this peak cannot be ruled out from the current observati ons. The h igh-energy component is more enhanced for the long- 
duration sources (Atoyan, Aharonian & Volk 1995; Kawan aka. loka fc Noiiri 2010). This is because the longer the duration 
of injection is, the larger fraction of fresh electrons/positrons are expected to reach the Earth without losing their energy so 
much during the propagation. MSPs are the continuously injecting sources with the duration as long as the cosmic age, so 
that the spectrum from them has nearly the same shape with the injection spectrum with the soft energy tail component. 
This is the difference from other sources such as young pulsars, whose typical duration is only ~ lO'^-lO^ yrs. Fitting result 
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of IVenter. Hajding &: Guillemot I2OO9') showed that the fraction of the pair starved MSPs is 25%, although they have only 



eight samples. Even if the fraction is 10%, the flux is ~ 20m~^ s~^ sr~^ GeV^ at 10 TeV. In this case, we can detect the 
electron/positron flux with near future missions such as CALET (we assume the geometrical factor times the observation 
time ~ 5 yrs as ~ 220 m^ sr days) because the predicted electron/positron flux is sufficiently large. It was considered that 
the number of astrophysical sources contributing to the above several TeV e nergy range is quite small according to the birth 



rate of the supernovae and t he canonical pulsars in the vicinity of the Earth (jKobavashi et al.ll2004l : iKawanaka, loka fc Noiiri 



2OI0I : iKawanaka et ahlboill ). However, we find that it is possible for multiple pair starved MSPs to contribute to the 10 TeV 



energy range in the electron/positron spectrum. Therefore, if the anisotropy of the observed electrons/positrons are weak in 
the 10 TeV energy range, we suggest that the pair starved MSPs may contribute to the spectrum significantly. 

Next, we also investigate the contribution of MSPs with high pair multiplicity to the observed cosmic ray electrons/ 
positrons. In this model, we assume that the injection function of these MSPs is power-law distribution with index a = 2, 
the cutoff energy £e,cut = ITeV and the minimum energy ee,min ~ IGeV. In this case, the pair multiplicity is 2000. Other 
parameters are the same values as in the case of the pair starved MSPs. We assume that the fraction of MSPs with high 
multiplicity is 100%. The results are shown in figure [5] for the electron/positron spectrum and in figure [3] for the positron 
fraction. Both figures show that electrons/positrons ejected from MSPs with high multiplicity partially contribute to the 



ical 


aulsars (Shen 197o|: Atovan. Aharonian & Volk 19951: Chi. Cheng & Youne 1996: 


""i->-i — „„„„ — „„„ 
Zhane & Cheng 2001: .Kobavashi et al. 


2004 


; Grimani 2007; Biisching et al, 20081: IProfumol 2012: Hoooer. Blasi & Seroico 


2OO9I: Yiiksel. Kistler & StanevI 2009: 


Malvshev. Cholis & Gelfand 20091: ICrasso et al.l 20091: Kawanaka. loka & Noiiril bOK 


: Hevl. Gill & Hernauistl 


20ld) would 



also contribute to the observed electron/positron spectrum. Note that even if other sources are dominant for the observed 
excess, the total spectrum added to the contribution of MSPs with high multiplicity does not significantly exceed the observed 
electron/positron spectrum. 



4 SUMMARY 

In this paper, we show the possibility that the cosmic ray electrons/positrons from MSPs would significantly contribute to 
the observed spectrum. Although MSPs have relatively low spin-down luminosity and low birth rate, the lifetime is so long 
that there are many active MSPs in the vicinity of the Earth. Furthermore, such a long lifetime source contirmously injects 
electrons/positrons after the nebula ceases expanding, so the adiabatic energy losses in a pulsar wind nebula region are 
negligible. The synchrotron cooling in the nebula is also small, so the high-energy electrons/positrons can escape the nebula 
without losing much energy. We calculate the diffusive propagation of high-energy electrons/positrons in the ISM taking into 
account the cooling via synchrotron emissions and inverse Compton scatterings, and predict their spectrum observed at the 
Earth. 

In the case of the MSPs with multiplicity that is unity, the typical energy of the electrons/positrons produced should be 
^ 50 TeV based on the assumption of the equipartition in the MSP wind region. Since the long duration of injection make 
a hard spectrum, the peak is enhanced at 10 - 50 TeV energy range. Even if the fraction of pair starved MSPs is as small 
as 10%, this peak would be detectable in the future missions such as CALET and CTA. Although a single young source can 
make the similar spectral feature in this energy range, in the case of pair starved MSPs the anisotropy of the electron/positron 
flux would be weaker because a number of sources contribute to it. If this peak is detected, that will be a great impact for on 
the studies of pulsar radio emission mechanisms because the existence of pair starved MSPs suggests that the radio emission 
mechanisms should be insensitive to the particle number de nsity down to sub-GJ num ber density. The detection also suggests 
that the current outer gap model should be modified because Wang fc Hirotariil (2011) suggested that most MSPs locate above 



the pair death line of the outer gap and the multiplicity is larger than unity. 

We also calculate the electron/positron spectrum from MSPs with high pair multiplicity. We suggest that if multiplicity 
of these MSPs is the order of ~ 10'', electrons/positrons from them partially contribute to the observed excess of the total 
spectrum and the positron fraction. 
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Figure 1. Electron/positron spectrum predicted from MSPs with the fraction of pair starved MSPs 100% (thin soUd Hne), and its sum 
(thick solid line) with the background (thin dashed line). The injection distribution function is the mono-energetic distribution ea. l|25|l 
with the energy £e = 50 TeV and the multiphcity k = 1. Data points correspond to measurements of ATIC (purple boxes, Chang et al. 
2008), HESS (light-green triangles and black triangles, Aharonian et al. 2008; 2009), PPB-BETS (yellow triangles, Torii et al. 2008b), 
Fermi (blue circles, Ackermann et al. 2010) and GRAPES (black circle, Kistler & Yiiksel 2009). We also show the total spectra from pair 
starved MSPs with the different fractions: 25% (thick dashed line) and 10% (dot-dashed line). We assume that the lifetime t = 5 X 10^'^ 
yr, the total energy -Erot = lO^'^ erg, the local birth rate R = 3 X 10~® yr~^ kpc~'^ and the fraction of the energy loss is 30%. 
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Figure 2. Electron/positron spectrum predicted from MSPs which have the power-law injection function with the cutoff energy £e,cut = 1 
TeV, the minimum energy £e,min — 1 GeV and index ct — 2 (thin solid line). The thick solid line shows the total spectrum. Por other 
parameters, we use the same values as in the case of Figure 1. 
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Figure 3. Total positron fraction resulting from the spectrum (solid line) that have the same parameters as in Figure 2, and the 
background (dotted line), compared with the PAMELA data as red circles ll Adrian! et al.l [20091 ). Note that the solar modulation is 
important below ~ 10 GeV. 
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